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REPORT 995 


BLOCKAGE CORRECTIONS FOR THREE-DIMENSIONAL-FLOW CLOSED-THROAT WIND 
TUNNELS, WITH CONSIDERATION OF THE EFFECT OF COMPRESSIBILITY 

By John G. Herriot 


SUMMARY 

Theoretical blockage corrections are presented for a body of 
revolution and. for a three-dimensional unswept wing in a cir- 
cular or rectangular wind tunnel. The theory takes account of 
the effects of the wake and of the compressibility of the fluid, 
and is based, on the assumption that the dimensions of the model 
are small in comparison with those of the tunnel throat. For- 
mulas are given for correcting a, number of the quantities, such 
as dynamic pressure and Mach number, measured in wind- 
tunnel tests. The report presents a summary and, unification 
of the existing literature on the subject- 

INTRODUCTION 

When a model is placed in a closed-throat wind tunnel 
there is an effective constriction or blockage of the flow at 
the throat of the tunnel. The effect of (his blockage is to 
increase the velocity of the fluid flowing past, the model; 
if the model is not too large relative; to the tunnel throat, this 
velocity increment is approximately the same at all points 
of the model so that the model is effectively working in a 
uniform stream of fluid the velocity of which is, however, 
greater than the frec-stream velocity observed at some 
distance upstream of the model. It is therefore necessary 
to correct the observed velocity, dynamic pressure, Mach 
number, and other measured quantities for the effect of 
this constriction. This correction is frequently called the 
correction for “solid blockage.” In addition to this solid- 
blockage correction, a correction for “wake blockage” is also 
necessary if the true dynamic pressure and Mach number at 
the model are to be determined. This wake blockage 
arises because the fluid is slowed down in the wake and con- 
sequently must be speeded up outside the wake. A further 
effect of the wake is to produce a pressure gradient which 
must be considered in correcting the drag coefficient. 

Formulas for the solid-blockage correction for a model 
mounted in a two-dimensional-flow wind tunnel are given in 
references 1, 2, 3, and 4. The first-order effects of the com- 
pressibility of the fluid on these corrections arc given in 
reference 5, as well as in references 3 and 4. References 
2, 3, and 4 consider also the wake-blockage correction, and 
reference 3 also considers the drag correction due to the 
pressure gradient caused by the wake. The formulas given 
in reference 3 for the solid- and wake-blockage corrections 
will usually be found most convenient whenever the engineer 
is confronted by a practical problem of determining the 
corrections for any configuration met in his experimental 
work. 


The solid-blockage correction for a model mounted in a 
three-dimensional-flow wind tunnel has been given in a 
number of different forms by different, authors. Not only 
do different authors give the correction for the same con- 
figuration in different forms but no one author gives formulas 
which are applicable to both fuselages and wings in tunnels 
of various shapes; for this reason, the engineer confronted 
with a correction problem may have to refer to several reports 
to get the complete solution of his problem. Moreover the 
modifications of the formulas for the first-order effects of 
fluid compressibility arc given incorrectly in some cases. 
References 1 and 2 give a formula for the solid-blockage 
correction for a body of revolution in a circular or rectangular 
tunnel for the case of incompressible flow'. In references 5 
ami 6 the effect of the compressibility of the fluid on this 
correction is discussed, but the result given is incorrect. 
Reference 4 gives a formula for the solid-blockage correction 
for a body of revolution and for a three-dimensional wing in 
a 7- by 10-foot wind tunnel. The modification for com- 
pressibility is correct for the wing but WTong for the body 
of revolution. Reference 7 gives a formula for the solid- 
blockage correction for a body of revolution in a circular 
tunnel, correctly taking account of the effect of the com- 
pressibility of the fluid. References 8, 9, and 10 give a 
formula for the solid-blockage correction for any body in a 
circular wind tunnel together with the appropriate constants 
for a body of revolution and for a rectangular wing having 
various span-to-diameter ratios; the modification of this 
formula to take account of compressibility is correctly given. 

It is clear that no one report gives all the necessary 
formulas together with the appropriate constants and com- 
pressibility modifications to enable the engineer to calculate 
the solid-blockage correction for any case with which he 
may be confronted. Moreover, when the results of two or 
more reports overlap, the forms arc frequently different so 
that it is not obvious whether the results are in agreement. 
It is the purpose of the present report to summarize and 
extend the results of the previously mentioned reports. 
Formulas are given for the calculation of the solid-blockage 
correction for a body of revolution or a three-dimensional 
unswept wdng in a circular or rectangular tunnel. These 
formulas contain t wo constants, one depending on the shape 
of the body and the other on the shape of the tunnel and 
the ratio of wing span to tunnel breadth. (This ratio may 
be taken to be zero for a body of revolution.) Values of the 
first constant for various bodies of revolution and for a 
number of frequently encountered wing-profile sections are 
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given. Values of the second constant for a circular tunnel 
and for rectangular tunnels of a number of commonly en- 
countered bread th-to-height ratios are given for various 
wing-span-to-tunncl-brcadth ratios. Some of these values 
have been taken from the previously mentioned reports; 
whereas others appear for the first time in the present, report. 
The discussion is limited to bodies centrally located in the 
wind tunnel. 

The wake-blockage correction for a model mounted in a 
three-dimensional-flow tunnel is given in references 4, 8, 9, 
and 10. For the case of incompressible flow the formulas 
are in agreement, but the modification to take account of 
compressibility given in reference 4 differs from that given 
in references 8, 9, ami 10. This matter is discussed in the 
present report. In reference 1 1 the corrections for the 
pressure gradient duo. to the wake arc given with the correct 
compressibility factors. 

All the final correction formulas together with directions 
for their use are given in the final section entitled “Conclud- 
ing Remarks.” Mathematical symbols are defined as intro- 
duced in the text. For reference, a list of the more important 
symbols and their definitions is given in appendix B. 

SOLID BLOCKAGE IN INCOMPRESSIBLE FLOW 

In studying the flow over a thin airfoil of small camber at 
a small angle of attack it, has been shown that the effects of 
camber and thickness may be considered independently. In 
treating the problem of wall interference, it is again con- 
venient to consider the thickness and camber effects sep- 
arately. The camber effect, as pointed out in reference 3, 
contributes nothing to the blockage correction. Conse- 
quently it suffices to determine the blockage correction for 
symmetrical bodies at zero angle of attack. This means 
that for wings it is necessary to consider only the base profile 
of the airfoil, the base profile being defined as the profile the 
airfoil would have if the camber were removed and the re- 
sulting symmetrical airfoil placed at zero angle of attack; 
bodies of revolution need be considered onty at zero angle of 
attack. 

RECTANGULAR TUNNEL 

Three-dimensional wing. — The blockage correction for a 
three-dimensional wing in a rectangular tunnel is considered 
in reference 4. Numerical values arc given only for a wing 
of 6-foot span in a. 7- by 10-foot wind tunnel. The method 
may, however, be applied to any rectangular tunnel and any 
span-to-breadth ratio. 

Consider a rectangular tunnel of height II and breadth B. 
-Suppose that the wing span 2s- is in the B-direetion so that 
2s/B is the span-to-breadth ratio. As in reference 4, let 
the wing be represented by a series of finite lines of sources 
and sinks. (See fig. 1.) The strengths of these sources 
and sinks are assumed to depend on the airfoil profile and 
the determination of these strengths, which is a two- 
dimensional problem, is -explained presently. If each line 
retains the same strength from one end to the other, the 
wing section cannot be exactly constant. It will thin 
down at the extreme tip and the plan form will not he 
exactly rectangular, but neither of these features is such as 
to detract from its usefulness for the present purpose, which 



Fihure 1 .— Image system for three-dimensional wing in rectangular tunnel. 


is to represent an actual wing sufficiently well to enable a 
calculation to be made of the velocity along the tunnel 
axis “induced” by images of the wing. It is this velocity 
induced hv the images which represents the effect of tin 
tunnel walls on the velocity at the model, and which is the 
solid-blockage correction. 

Consider the image line source CD of strength Q per unit 
length. The velocity potential of a three-dimensional 
source of strength Qby (volume per unit time) is — Qby/Awr. 
It follows that the component velocity along the tunnel 
center line induced at A by the source element Qby is 


bu A 


Q(J 

47 T /' 3 


Dj 


where y and r arc defined in figure 1. Putting r 2 = 
(mB-y) 2 +n 2 H 2 -\- g 2 and integrating from — s to s gives 


Aw, i 


Qg I" mB+s 

4 7r {n 2 H 2 + g 2 ) _ y n 2 IP + g 2 -\- ( mB + s) 2 


mB—s 

yTm 2 +g 2 +(mli-C 




(1) 


for the single line source CD. Now CD is one image of one 
finite line source used to represent the wing in the tunnel. 
In order to find the velocity induced by all the images of this 
particular line source, it is necessary to add the results 
obtained from equation (1) by giving m and n all positive 
and negative integral values except m — n= 0. Thus for a 
single line source 

.. v” 1 T mB+s 

1 4 7r n 2 H 2 -yg 2 \_ A ]n 2 H 2 -Yg 2 -\-{mB-)rC) 2 


mB—s 

^W+g^+JmB-s) 



( 2 ) 
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where the prime denotes that the. term ni—n = 0 is to be 
omitted from the summation. Equation (2) may be re- 
written 

.. _ Q'J# ( (J * /i\ 

1 2 ir IP a \/J ! B' II ) 

where 
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Now a and r depend only very slightly on ///// and so it 
usually suffices to take {llll= 0 in the evaluation of these 
quantities. Any wing profile can be represented by a suitable 
distribution of sources and sinks along the chord. It follows 
that, the. total induced velocity due t,o I, he wing images is 
obtained by summing over this distribution and is approxi- 
mately 


A ,U' 


JZQJL 

(lil-iyi* 



* IT) 
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It may lie noted here that setting <///-/=() in the evaluation 
of a and r is equivalent to representing the wing by a line 
doublet of strength "ZQy (analogous to references 1 and 2) 
instead of by a distribution of sources and sinks. 

The quantity 2 Qy of equation (5) can be determined 
approximately from the wing profile and this determination 
is a two-dimensional problem. From reference 2, but with 
the notation of reference 3, there is obtained 

ZQ(/=^Ac 2 IJ’ (6) 

where 

c airfoil chord 

U' apparent frcc-stream velocity at airfoil as determined 
from measurements taken at a point far ahead of 
model 

A a factor dependent on shape of base profile 

Substitution of equation (G) into equation (5) yields 

A JP 2 $c* P ' 2 ( « B\ 

U' (Bliy 12 16 \ ’ B’ II J 


2 set A / s B\ 

C vz 16 tie T Y ! B’h) 


(7) 


where 

t maximum thickness of airfoil 
(J cross sectional area of tunnel 

If V denotes the volume of the wing, then V’=2sc/k 1 where 
ki < 1 depends on the shape of the base profile. Equation (7) 
may be rewritten 


or 


where 


I t is clear that r depends only on the tunnel shape and the 
wing-span-to-tunnel-brendth ratio; whereas K l: K 2 depend 
only on the shape of the base profile. 

The factor A can be determined for any base profile from 
the relation (references 2 and 3) 

A =~ j’' V ~ /l -E Y+(dy,/dxyd (f ) ( 1 3) 

where 

y, ordinate of base, profile at chordwise station x 
di/t/dx slope of surface of base profile at x 
P base-profile pressure coefficient at x in an incompres- 
sible, flow 

Values of A for a number of base profiles are given in refer- 
ence 3. Thus the value of K? can be calculated from equa- 
tion (12). The value of K t is immediately found from the 
area of the base profile which may be calculated, for example, 
by a numerical integration from the ordinates of the base 
profile. As soon as *ci is known, K, can then be calculated 
from equation (1 1). The evaluation of r requires the sum- 
mation of the infinite series in equation (3). The summation 
of this series is explained in Appendix A. 

Values of r for rectangular tunnels of various breudth-to- 
height. ra tios and for various wing-spa n-to-tunnel-breadth 
ratios are given in table I and figure 2. Values of K\ and K 2 
for various base profiles are given in tables II and III. The 
choice between the two formulas (8) and (9) is entirely a 
matter of convenience and should be decided in the light of 
the available da ta. 

Body of revolution. — The blockage correction for a body 
of revolution in a rectangular tunnel is considered in refer- 
ences 1, 2, and 4. In reference 4 numerical values are given 
for some average streamline body of revolution, in a 7- by' 
10-foot wind tunnel; whereas in references 1 and 2 numerical 
values are given for prolate spheroids and Rankine Ovoids 
in square and duplex (breadth equal to twice the height) 
wind tunnels. Either the method of reference 4 in which 
the body' of revolution is represented by' a suitable distribu- 
tion of sources and sinks along its chord or the method of 


A ,(/' Ki t V 
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A ,U' K 2 t2scI 
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TABLE I.— VALUES OF r FOR VARIOUS TUNNEL SHAPES 
AND CONFIGURATIONS 




llwlv of 

3-dimensional wing with span-to-breadth ratio 

Timnrl shai>e 


revolu- 

tion 

fr n 

rr 0 ™ 

f=0.60 

1=0,6 

?->» 

Circular 

Square 


0.797 

.812 

0. 797 
.812 

0.812 

.SIS 

0. 828 
. 886 

0. 859 
.874 

0. 951 

mi i=ioi7 

. m 

. 868 

. 864 

.866 

.884 

.915 


mu — ha 

.945 

. 946 

.911 

. 930 

. 923 

. 937 


mn=2 

1.028 

1.028 

1.017 

.990 

. 957 

. 952 

Rectangular... 

mi i=H2 

1. 729 

1.729 

1.680 

1 . 430 

1. 204 

1. 150 

m n=2i7 

1.729 

1. 729 

1 . 783 

1.806 

2. 195 

2. 555 


Ilf 1 1=1/2 

1.028 

1.028 




1 . 434 


mi i=o/ to 

.922 

.922 

. 932 

.977 

1.063 

1.214(1 

, 

!«///= 7/10 

.868 

.863 




1.110 


TABLE II.— VALUES OF A', FOR VARIOUS BASE PROFILES 


\ 








\ 



Con ven- 

NACA low-drag sections 

\ 


Joukow- 

t.ional 





, \ 
r o.xx\ 

Ellipse 

ski 

section 

NACA 

sections 

ooxx 

15-0XX 

64-OX X 

65-OX X 

55-OX X 

o. or. 

0. 938 


0. 941 

0. 938 

0. 952 
. 901 

0.955 

0.987 

.09 

. 955 

0. 991 

. 972 

. 952 

. 907 

. 984 

. 12 

. 993 

1.016 

1.005 

. 997 

. 987 

. 989 

1.000 

.15 

1.019 

1.045 

1 . 035 

1 . 028 

1.019 

1.020 

1.0.32 

. 18 

1.040 

1.058 

1 . (Mil 

. 1.052 

1.047 

1.051 

1.057 

.21 

1.072' 

1.083 

1 . 090 

1.090 

1.073 

1 . 079 

1.079 

.25 

1. 108 


1. 128 

1. 128 




.30 

1 . 152 






. 35 

1.195 







.50 

1.329 







1.00 

1.772 














TABLE III.— VALUES OF K 2 FOR. VARIOUS BASE PROFILES 


\ 












Conven- 

NACA low-drag sections 

X 

Ran- 


Joukow- 

tionnl 





X 

kinc 

Ellipse 

ski sec- 

NACA 





i X 

=o.XX\ 
c \ 

Oval 


tion 

sections 

OOXX 

IMXX 

54-OX X 

65-OX X 

65-0X X 

0. 00 


0. 737 


0.044 

0. 690 

0.015 

0. 632 

0. 679 

.09 

0.913 

.758 

0.599 

. 665 

.708 

.011 

. 030 

. 673 

. 12 

.928 

.780 

.015 

.587 

.734 

. 024 

.041 

.084 

. 15 

. 935 

.800 

.033 

.70S 

.750 . 

.540 

. 057 

. 098 

. 18 

. 953 

.822 

. 052 

.727 

.781 

. 054 

.073 

. 712 

.21 

.901 

.842 

. 070 

.740 

.802 

. 005 

.686 

.723 

.25 

.979 

.870 

. 592 

.771 

.830 

.678 

. 097 

.732 

.30 

1.002 

. SKI5 

. 725 






.35 

1.043 

. 940 

.763 






. 50 

1. 170 

1.044 

.876 






1.00 

1.392 

1.392 








references 1 and 2 in which the body is represented by a 
doublet of suitable strength at its center may be used to 
obtain results for any streamline body of revolution in any 
rectangular tunnel. Since the method of reference 4 was 
used for the three-dimensional wing, it is instructive to use 
the doublet method of references 1 and 2 for the body of 
revolution case, although both methods give the same results. 

Consider a body of revolution of maximum thickness t and 
length c centrally located in a rectangular tunnel of breadth 
B and height II. As in reference 2 the body may be repre- 
sented by a doublet of strength b given by the equation 
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(n) Vnriation of r with rat io of wing span to tunnel breadth, 2 »//?, for various tunnel sha|>cs. 
Figure 2.— Values of r for various tunnel shapes and configurations. 
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(b) Variation of r with ratio of rectangular tunnel breadth to height, BUI, for a body of 
revolution and for wings of various span-to-tunnel-brcadth ratios. 


Figure 2.— Concluded. 


the summation being taken over all positive and negative 
integral values of m and n except m=n= 0. If g/H=0 
equation (3) may be rewrit ten 


<7 


/ s B\ B 1 r(rn+s/B)F mn —(rn—slB)E mn 
V’B’hJ 2*^n 2 L E mn F mn 


= £' 


2 m 

Em n E m n [(»l + S/B) F m n + ~ S/B ) E m „] 

(16) 


where the quantities E mn and F mn , which arc introduced for 
convenience, arc defined by the equations 


E mn = V n 2 +{m + sjB) 2 (B/H ) 2 



(14) 


F m „ = V»i 2 + (m - s/mB/I-D 2 

It follows that 


where A is a constant depending only on the shape and fine- 
ness ratio of the body. The velocity induced at the model 
by the tunnel walls is the same as that induced by a doubly 
infinite array of images of the doublet and is given by 


<7 (<>, Or S' [ n 2 +m 2(BIH )*\ 3/2 

Substitution of equation (17) into equation (15) yields 
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If is replaced by its value from equation (14) and equation 
(4) is used, there is obtained 


A ,(/' ir V2 A/, 3 

U' 8 (I3II) 312 T 



~C 3 ’ 2 S c T 



If I' denotes the volume of the body of revolution, then 
I'=K 2 ct 2 where k 2 <it / 4 depends on the shape of tins meridian 
section of the body. (For a right, circular cylinder whose 
meridian section is clearly a rectangle k 2 =7t/4.) Equation 
(18) may be rewritten 



A JJ’ K 3 tV 
LT C 312 

(19) 

or 

A ,(/' K t rct 2 
U' C 3 ' 2 

(20) 

where 

° 4) 

(21) 


7 T 3,2 \t. 1 

8 c « 2 

(22) 


w 3/2 \t 

A)= T7 

(23) 


TABLE IV.— VALUES OF K, FOlt VARIOUS BODIES OF 
REVOLUTION 


\ 

i \. 

* \ 

Prolate 

spheroid 

Rankine 

Ovoid 

NACA 

111 

NACA 

133 

Fuhr- 

mann 

source- 

sink 

body 

Sym- 

metric 

Fuhr- 

mann. 

source- 

sink 

body 

(3, 6. 
-1,0) 
body of 
reference 
14 

0. 08 

0.900 

0. 913 

0. 902 





. 10 

.905 

.923 

.909 


BHtHh 



. 12 

. 910 

.932 

. 910 

B RwB 

bhhbh 



. 14 

.917 

. 94 1 

.924 





. 10 

. 924 

.949 

.931 





. 18 

.931 

.957 

.939 


0. 933 

0. 932 

0. 925 

. 20 

. 938 

. 904 

.948 





.24 

.954 

.979 

.904 





.30 

.980 

1.000 

i • 

0. 983 




.40 

1.025 


V ' ' 



■ 



1. 072 








1.329 








TABLE V.— VALUES OF K, FOR VARIOUS BODIES OF 
REVOLUTION 


\ 

t \ 

’ \ 

Prolate 

spheroid 

Rankine 

Ovoid 

NACA 

111 

NACA 

133 

Fuhr- 

mann 

source- 

sink 

body 

Sym- 

metric 

Fuhr- 

mann 

source- 

sink 

body 

(3. 0, 
-1.0) 
body of 
refer- 
ence 14 

0.08 

0. 472 

0. 070 

0. 405 





. 10 

.474 

. 008 

. 409 





. 12 

.477 

. 005 

.413 





. 14 

.480 

. 002 

.417 





. 10 

.484 

. 059 

. 422 





. 18 

. 488 

. 055 

.428. 


0.417 

0.414 

0. 430 

.20 

.491 

. 050 

. 433 





. 24 

. 500 

. OH) 

. 445 





.30 

.514 

. 024 






.40 

.537 

. 007 






.50 

.502 

.009 






1.00 

• 69C 

.690 















It should be noted that r for the case of the body of revo- 
lution is the same as for the limiting case of a wing when the, 
span approaches zero, As pointed out- in reference d, A may 
he calculated for any body of revolution whose pressure dis- 
tribution is known. The necessary formula is 

(j)'P. (fTvr^Vi +<<»<.)’<< (I) wo 


where 

y radius of body at chordwise station x 
dyjdx slope of meridian section at a: 

P pressure coefficient at x in incompressible flow 
Values of X for prolate spheroids and Rankine Ovoids are 
given in references 1 and 2. As soon as X is known for a 
body, K, can he calculated at once from equation (23). . The 
value of *2 may be found from the volume of the body of 
revolution which may be calculated for example by a numcr- 



NACA III ; t/c -0.20 



Prolate spheroid ; t/c - 0.20 





h/AC A III ; t/c -0.18 



Fuhrmonn source-sink body ; t/c =0.18 



Symmetric Fuhrmann source- sink body, t/c *0,18 



(3,6,-!,OJ Body ot reference I4\ t/c =0.18 
Figure 3.— Sample bodies of revolution. 
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ical integration, ami K 3 can be calculated from equation (22) 
as soon as k 2 is known. 

Values of t for rectangular tunnels of various brendth-to- 
lieight ratios are given in tabic L and figure 2. Values of 
Ii 3 and /v'( for various bodies of revolution arc given in 
tables IV and V. Some of these bodies are drawn in figure 
3. Unfortunately these tables arc rather incomplete and 
moreover fuselage shapes used in practice are extremely 
varied. However, in table IV it is observed that the values 
of K 3 do not de])cnd very strongly on the. shape of the body, 
although they do depend on the. thickness ratio. For this 
reason it appeal’s that for most, fuselages it. will be sufficiently 
accurate to use the values of K 3 given for the NACA .'Ill- 
bodies. As the values of K t given in table V are more de- 
pendent on the body shape, it is recommended that equation 
(19) and table IV be used in preference to equation (20) and 
table V whenever possible. 

CIRCUI.AH TUNNEL 

Body of revolution. — It is convenient to consider the case 
of a body of revolution before considering the case of the 
three-dimensional wing because, more attention has been 
given to the former by other authors. (See references 1, 2, 
7, 8, 9, and 10.) It will now he. shown that the blockage 
correction is again given by equations (19) and (20) whore 
t has a value appropriate to a circular tunnel. Since K 3 and 
K t depend on the model and not bn the tunnel, they are 
still- given by equations (22) and (20). 

The most convenient, starting point is the- formula of 
references 1 and 2, namely, 




3/2 


(25) 


where S m is the maximum cross-sectional area of the model. 
It is only necessary to note that 


where 

factor depending on model shape 
t v factor depending on tunnel shape 
V volume of model 
B diameter of wind tunnel 

It is of interest to show that, equations (19) and (20) can la 
deduced also from equation (27) showing the latter to la 
equivalent to equation (25). From reference 8, there i> 
obtained 


X _liL x 

Xr ~4 U X 


Ty=- T 

ir 


Since also C=v/F/4 equation (27) yields at once 


-AiLT 

U 


U' /it t 3 \ / 4 \ „/V /2 1 \ Xf 

-“(i V V(i V V \T &*)=’ 7 T 


™ Cl 2 

C w 


This is (he same as equation (20) which yields equation!' 
(19) and (20) directly. 

Three-dimensional wing. — The blockage correction for a 
three-dimensional wing in a circular tunnel is given in refer- 
ences 8, 9, and 10, the formula being of the same form as for 
a body of revolution in a circular tunnel, namely, equation 
(27) where r v and \ v have values appropriate to the three- 
dimensional wing. From reference. 8 (using the. notation of 
reference .3 instead of reference 2) there is obtained 


x - 7r \ 

A V jj ~~y~ A 


(28) 


4 

T|’=- T 
7 T 


(29) 


Since also C—B^/d equation (27) yields at once 


S m .= 7rt 2 /4 

1 = 

It follows at once that 

A, U' , t F /2 I 3 \t tt^ 2 c / 2 

w rX S C 312 7 c S C Vi 


( 20 ) 


Substitution from equations (22) and (23) reduces this to 
equations (19) and (20), respectively, ft should' he noted 
that r of references t and 2 is identical with r of t lie present 
report. 

The value of r is given in table 1 and figure 2. Values of 
K 3 and K t are given in tables IV and V. Again equation 
(19) is preferable to equation (20). 

The result of reference 7 fails to take into consideration 
the body 7 shape and so it is not very useful except for less 
exact calculations. The Jesuit of references 8, 9, and 10 is 
presented in a different- form, namely. 


A K' 
U’ 


— X v T v 


w 


A|(/' {v 2xc 2 ,V'4 \ ir /V' 2 1 \ 2 set A r V2 

U' ~VS l’ V \ 8 C 3 ' 2 ) C 3 ' 2 tfc 10 r 

This is the same as equation (7) which leads directly to 
equations (8) anil (9) with the same definitions of K t and 
Ii, given in equations (1 1) and (12); r is, however, given In- 
equation (29) where 7>is obtained from refei’onces 8, 9, or 10. 

d ims equations (8) and (9) may also he used for circular 
tunnels provided only that the appropriate, values of r arc 
usedi Values of r are given in table I and figure 2. Values 
of K\ and K 2 are given in tallies II and 111. 

SOLID BLOCKAGE IN COMPRESSIBLE FLOW 

In the preceding section the solid-blockage corrections 
have been determined under the assumption that the fluid 
is incompressible. It is now necessary to determine the 
modifications required in these formulas to take account of 
•the ed’oets of the compressibility of the fluid. Tbe methods 
of refeiences 12 and 13 arc very convenient for this purpose. 
As the required modifications arc given incorrectly in refer- 
ences 5 and 6, partially incorrectly in reference 4, and 
correctly in references 7, 8, 9, and 10, it appears worth while 
to give some discussion of the matter; 


(27) 
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loir (lie purpose of deducing the properties of a compressi- 
ble How from (host 1 of a corresponding incompressible (low 
the so-called “Extension of t.lic Prandt.l Rule,” which was 
first, given in reference 12 and repeated as Method I V in 
reference Id, is prohahly of most general application; hut 
the other methods of reference Id are sometimes more con- 
venient for certain problems. The Extension of the Prandt.l 
Rule may he expressed in the following manner: 

The streamline pattern of a compressible (low to he cal- 
culated can he compared with the streamline pattern of an 
incompressible flow which results from the contraction of the. 
y and 2 axes including the profile contour by the factor 
-\/l— A/ 2 (iV/=frcc-sfream Mach number) (x axis in the 
direction of the free stream). In the compressible (low the 
pressure coefficient as well as the, increase in the longitudinal 
velocity are greater in the. ratio 1/(1 —AT 2 ) and the stream- 
line slopes greater in the ratio 1/yl— At 2 than those at- the 
corresponding points of the. equivalent incompressible flow. 

Since formulas (S) and (19) for the three-dimensional wing 
and the body of revolution have the same form and also 
apply to both rectangular and circular tunnels, it suffices to 
determine the modification due to compressibility for them. 
Let the subscript c refer to compressible flow and the sub- 
script i to the corresponding incompressible How. If V c is 
tlu! volume, of the model in the compressible How, then T' ( , 
the volume of the model in the corresponding incompressible 
How, is given by 

v t =[i-(Aiy}V c 

where AT is the apparent frec-stream Mach number at the 
model as determined from measurements taken at a point 
far ahead of the model. Also if C c is the cross-sectional 
area of the tunnel in the compressible flow, then Ct the 
cross-sectional area of the tunnel in the incompressible flow 
is given bv 

t is unaffected by the transformation and the effect on A"j 
and K 3 is sufficiently small that it may be neglected. From 
the Extension of the Brand tl Rule it follows that 


(*V'\ _ 

1 , 

(AJJ'A 

V u' )c 

■ 1 -(Al'f 1 

V U' )i 


1 K)TVi 1 KjT V c 

l -{ Al'Y c\ m [1 — ( i !/') L ’] :,/2 C 3 ' 2 

where j denotes the numbers 1 or 3. Since equations (S) 
and (9) are equivalent as are also equations (19) and (20), 
it follows that in all cases it is only necessary to multiply 
the blockage corrections given by these formulas by 
[1 —(M') 2 ]~ 312 in order to take account of the compressibility 
of the fluid. 

As a check it is useful to determine the compressibility 
modification for the case of a body of revolution in a circular 
tunnel by Method II of reference 13, since the derivation 
is so simple by this method. Both the body shape and 
the longitudinal velocities arc the same in the corresponding 
compressible and incompressible (lows; only the tunnel 
dimensions are altered by the factor -y 1 — (.'V/') 2 so that 
— (At’) 2 ]C c . There is obtained 


r 


A't T Vi 1 K 3 t V c 

C 312 [\-(M') 2 1 3/ * Q, 3n 


as before. 


WAKE BLOCKAGE IN COMPRESSIBLE FLOW 

The blockage due to the wake of a model in a two-dimen- 
sional-flow tunnel is discussed in some detail in reference 3. 
Much of this discussion is applicable without change to the 
case of a, three-dimensional (low tunnel. The. fundamental 
idea of replacing the model and its wake by a. source (in this 
case a three-dimensional source) of suitable strength located 
at the position of the model can be used again and in fact the 
determination of the source strength Q can be carried out in 
exactly the same manner. The result is identical with that 
of reference 3, namely, 


Q- 


p'U'Cd'S 


[H-(t-1)(M' 


(30) 


where 

Q mass flow of source rather than volume flow as used 
previously 

P mass density of fluid at point far upstream 
G a uncorrected drag coefficient referred to apparent dy- 
namic pressure q' 

S area on which drag coefficient is based 
7 ratio of specific heat of gas at constant pressure to 

specific heat at constant volume (~f^ 

In equation (30) powers of At' higher than (At') 1 have been 
neglected. 

Consider now a reetagular tunnel of height It and breadth 
B. The tunnel walls are replaced by a doubly infinite array 
of sources of strength Q at, distances inB to the side and nil 
above and below the position of the model. 

By making use of Method II of reference 13, it, is readily 
shown that a three-dimensional source of strength Q (mass 
per unit time) in a uniform flow of compressible fluid will 
induce at the point the coordinates of which are x, y, z rela- 
tive Lo the source a stream wise velocity. 


Au — 


Q. 


4irp [x 2 -f-(l — M t )(y 2 -\-z 2 )] 312 

where the uniform flow is in the x direction and p and M are 
the density and Mach number, respectively, of the undis- 
turbed stream. 

It follows that the streamwise velocity A 2 U' induced at a 
point on the center line of the tunnel by the entire system 
of images is 

<3 


AoU'-- 


4 7 rp 


7 S' 


{z 2 + [ 1 — (At') 2 } (m 2 JP + n 2 !! 2 )} 


3/2 


where p' and At' are the density and Mach number of the 
undisturbed How in the tunnel and the summation is taken 
over all positive and negative integral values of m and n 
except ?« = « = 0. The velocity induced by the image 
sources at an infinite distance upstream is 

Q 


(A 2 t/')- 


lim _ y >' 


yKm 2 A 2 -H/3//q}w 



778 


REPORT 995 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


But. conditions far upstream must remain unchanged and to 
achieve this it is necessary to counterbalance this velocity 
by superposition of a uniform (low of equal magnitude but 
opposite sign. The addition of this flow at all points in the 
field will result in a speeding up of the general flow at the 
position of the airfoil by the amount 


Aa U'=lim -r-^-7 S' 

4 jrp 


{x 2 + \\-{M ')'i (m*JP + n 2 !! 2 )} 3/2 


(31) 


The summation of this series can he obtained by the following 
artifice. Substitution of equation (30) into equation (31) 
and setting- M' = 0 gives 


ft ' o / / ' 

(Ajj'U^hni XT 


Sir ^ {x 2 +m 2 B 2 +n 2 H 2 ) 312 
But according to references 4, 8, 9, and 10 
’'A,(/'\ 1 CVS 


Maf/'X _1 
V U' ), 4 


BI1 


Comparison of equations (32) and (33) shows that 

X 2 7T 


lim S' 7 


(x 2 + m?B 2 +n 2 HY 2 BH 


( 32 ) 


(33) 


(34) 


WAKE PRESSURE GRADIENT IN COMPRESSIBLE FLOW 


The effect of the pressure gradient caused by the wake is 
discussed in reference. 3 for compressible flow in a two- 
dimensional wind tunnel and in reference 1 1 for compressible 
flow in a circular wind tunnel. The method of reference 11 
is equally applicable to the case of a rectangular wind tunnel 
if the appropriate value of t (table. I) is used. 

The longitudinal velocity increment due to the effect of 
the tunnel boundaries on the source used to simulate the 
wake has an approximately 7 linear gradient in the stream 
direction at the model location. This linear gradient in the 
velocity 7 is equivalent to a linear gradient in the pressure as 
is easily seen from the approximate relation. 



A a £/' 

u' 


In reference 11 it is shown that the gradient in &?U' for a 
source in a wind tunnel is identically 7 equal todhe value of 
A 2 U' for a doublet in a wind tunnel. In the notation of 
the present report there is obtained 



For compressible flow this becomes 


If in. equation (34) B and II arc replaced by yd — {M') 2 B 
and VT —(M') 2 II, respectively, there is obtained 


</£__ ,. i+(7-i)(MT b Cj/S 
dx ~ q [1-(M') 2 ] 3/2 \ 4 T C 3 ' 2 


l !"l ^ {x 2 +[l-(M') 2 } (m 2 B 2 +n 2 H 2 )} 3,2 ~ [l —(M') 2 ]BH 

(35) 

Substitution of equations (35) and (30) into equation (31) 
yields 

A JU’ l+(y-l )(M') 2 C,/S \+0A(M') 2 C,>’S , ntn 
U' ~ 1 -(My 4BH 1 — {M') 2 4 C K ' ’ 

on setting 7=1.4. The preceding discussion is for a rectan- 
gular tunnel, but in references 8, 9, and 10 the same formula 
is given for a tunnel of any 7 shape for the case M' — 0. Con- 
sequently 7 , equation (36) may 7 be taken to hold generally 7 for 
a tunnel of any shape. 

In reference 4 the wake-blockage correction is given cor- 
rectly 7 for incompressible flow but the modification to take 
account of compressibility is given incorrectly, as it is deter- 
mined from reference 5 which, as pointed out in reference 13, 
is incorrect. The effect of compressibility 7 on the wake- 
blockage correction is determined in references 8, 9, and 10 
by means of the correct form of the Extension of the Prandtl 
Kule. but there is some question whether this rule is applica- 
ble to this problem: It appears that, in using this method 

to determine the effect of compressibility on the blockage 
correction, the effect of compressibility on the source strength 
Q given by. equation (30) is overlooked. This explains the 
discrepancy 7 between equation (36) and the formula of refer- 
ences 8, 9, and 10. 


, 1+0.4(M') 2 b Oo’S 
~ 2 [l-(M') 2 ] 3 ' 2 V 4 T C 3 ' 2 

The increase in the drag resulting from^this pressure 
gradient is equal to the product of the pressure gradient 
by the sum of the actual model volume and the virtual 
volume (reference 2). It should be noted that the constants 
A', and K 3 of tables II and IV are equal to V0r/4) times 
the ratio of the sum of the actual model volume and the 
virtual volume to the actual model volume, these quantities 
being calculated under the assumption of incompressible 
flow. Thus the increase in drag coefficient, caused by the 
pressure gradient is given by 


AGV = 

1+0.4(M') 2 K,tV w C d ' 
[1— (M') 2 ] 3/2 C 312 

(37) 

for the wing, and 



A ft > — 

1+0.4(M') 2 AVfiW 
[1— (M') 2 ] 3/2 C 3/2 

(38) 


for the body 7 of revolution. 

It is pointed out in reference 11 that the virtual volume is 
altered by the compressibility of the fluid. Thus equations 
(37) and (38) can be slightly 7 improved if K\ and K 3 appearing 
therein are corrected to take account of this effect. The 
necessary modification is made by 7 replacing A', and K% in 
these equations by K\ p and K 3p , respectively 7 , where 

K]p = V? [ 1 +h(M ') (\/^ ~ 1 )] ’ ^' = 1 ' 3 (39 > 
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Figure 4 .— Linear compressibility correction factor for virtual volume. From reference 11 . 

where h (M') is given by figure 4 which is reproduced from 
reference 11. It should he noted that, the improvement in 
modifying K\ and A : 3 for compressibility will be small es- 
pecially in the ease of Zv 3 for the body of revolution. It is 
important, only for quite high Macli numbers. This addi- 
tional refinement, was not made in reference .‘3. 

It appears that a similar compressibility modification of 
A, and A' 3 in the formulas for the solid blockage should be 


made, but the exact modification required is not known. 
However, it is behoved to be small. 

CORRECTION OF MEASURED QUANTITIES 

The true velocity U at a model consisting of a body of 
revolution and wing can be obtained from the apparent 
velocity U' bv applying the solid-blockage and wake- 
blockage corrections. The true velocity may be written in 
the form 

U — U'{\ +A) (40) 

where 

A = K w + Kt> + A+ (41) 


In equation (OS) K n is the solid-blockage correction due to 
the wing and is given by 

1 AjT V w 


A„= 


or 


A»= 


11 -{M')Y 2 C 3n 

.1 A2l’2.SCdtt 


[1— (M')T 


C m 


(42) 


(43) 


K b is the solid-blockage collection due to the body of 
revolution, being given by 

1 A 3 t V b 


A 


or 


A 


* [i -(Myf 2 c m 

1 IUTC b t b 2 


^ [1-(M')T /2 O 312 

K, ak is the wake-blockage correction and is given by 


AT, ft 


1 + 0. 4 (A/') 2 Co'S 

i -(My 4 c 


(44) 


(45) 


(46) 


It is evident that a correction to the apparent velocity in 
a compressible llow implies corrections also to the apparent 
density, dynamic pressure, Reynolds number, and Mach 


TABLE VI.— COMPRESSIBILITY FACTORS FOR CORRECTION EQUATIONS 
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number. These corrections are readily obtained on the 
basis of the usual assumption that the flow is adiabatic. It 
is assumed that the correction terms are small compared 
with unity, so that squares and products of these terms may 
be neglected. The analysis follows the lines of reference 3, 
and it is not necessary to repeat (lie details here. The fol- 


lowing equations arc obtained: 

P = P '[l~(M7/lj (47) 

4 = (/{l + [2-(i\/') 2 l/v} (48) 

/>’=/."{ 1 +[1 -0.7(A/') 2 ]A} (49) 

M= M'{ l +[ l +0.2(M') 2 ]/t} (50) 


The drag coefficient must be corrected for the effect, of the 
pressure gradient, due to the wake as well as to refer it to the 
correct dynamic pressure. There is thus obtained 

C D =Cn'{l-[2-(My]K}-AC D J-AC Db ', (51) 
where AC D J and A6 'd/ are given by equations (37) and (38). 

Numerical values of the functions of M' which appear in 
these equations are given in table VI. 

CONCLUDING REMARKS 

Data obtained from tests of three-dimensional models, 
which are small relative to the wind-tunnel dimensions, can 
be corrected for solid and wake blockage and for the pressure 
gradient due to the wake by means of the following equations: 


U = U'{\+K) 

(40) 

2 = «'{1-H2-(MW} 

(48) 

/f=/r{l-Hl-0.7(A'/') 2 ]A'} 

(49) 

M = A/Ml +[1 +0 .2 (M') 2 ]K} 

(50) 

C D = Co'{ 1 -[2-(A/') 2 )/y } -*C D J- 

A6V (51) 

In the preceding equations A' is obtained from the following: 

A — AM AM- K lr k 

where 

(41) 

r- 1 K\TV W 

u [ 1 — (A/') 2 ] 3/2 C 3 ' 2 

(42) 

or 


rr 1 AV^VCuhr 

(43) 

. [1 — (A/') 2 J 3/2 C 3/2 


1 A’ 3 t ly 

u- [l -(iV /') 2 ] 3/2 C 3/2 

or 

j- 1 A 4 T C„t 

71i >-[l_(j\/')*]3/2 (73/2 

1+0.4(A/') 2 C n ’S 
i -(My 4 c 

, l+0.4(A/') 2 /t,r V„C D ' 

U-(M')T /2 C 3/2 

, p , 1 -4- O^AA) 2 K 3 t \'t)Ci) 

[ i — (d/') 2 ] 3/2 • C 312 ~ 

In these equations r is a factor which for a body of revolu- 
tion depends only on the shape of the tunnel; whereas for a 
three-dimensional wing r depends on the ratio of the wing 
span to the tunnel breadth as well as on the tunnel shape. 
(The wing span is in the direction of the lmmcl breadth.) 
Values of r arc given in table I and figure 2. The constants 
Ki and K 2 for t lie three-dimensional wing depend only on 
the wing base profile shape and can be calculated by means 
of equations (11), (12). and (13). Values of A) and A' 2 for 
a number of wing profiles are given in tables II and 111. 
The constants A'j and 7v'< for the body of revolution depend 
only on the shape of the body and can be calculated by means 
of equations (22), (23), and (24). Values of these constants 
for a number of body shapes are given in tables IV and V; 
some of these shapes are drawn in figure 3. Since these 
tables are incomplete and fuselage forms are not standard- 
ized as arc wing sections, it is recommended that the values 
of Kt given in table IV for the NACA 111 series of shapes 
be used for any fuselage shape which does not differ too 
greatly from an NACA 111 shape. This implies that equa- 
tion (19) and table IV should be used in preference to equa- 
tion (20) and table V whenever possible. Numerical values 
of the functions of A/' which appear in the correction equa- 
tions are given in table VI. 

The constants K, and I\ 3 appearing in equations (37) and 
(38) may be modified for compressibility bj' means of figure 
4. The modified values of Ivi and K 3 are to be used in 
equations (37) and (38) onlv and not in equations (42) and 
(44). 


Ames Aeronautical Laroratorv, 

National Advisory Committee for Aeronautics, 
Moffet Field, Calif. 
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APPENDIX A 

SUMMATION OF THE INFINITE SERIES FOR <r(0, s IB, Bill) 


From equations (3) and (16) it is seen that. <r(0, s/B, B/H) 
may be written in the alternative forms 


■’H-D-s s ' K“lroTOr) < A » 


or 


^ (°' B’ h) Ji m n F m „ l(m + ,•/. 


•Jin 


IB) I 1 m „ T {in si B)E nn ] 
(A2) 


where the summation is taken for all positive and negative 
integral values of in and n except m—n = 0 and the (|unntit'ics 
E mn and F mn , which arc introduced for convenience, are 
defined by Lhe equations 

E mn = V« 2 -f-(m +s/B) 2 (B/H) 2 


I’m n = A 2 + (m - s/mBIH ) 2 

It is possible to sum the series for <r(0, s/B, B/H ) exactly 
onlj r when s/B=~. In this case equations (Al) and (A2) 
yield 




= 2 y) — li m 

n = \ fl m—*o 
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in — - 


+ 2 S 
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2 ( m+ |) 


J n i + (m+ Y y C BUT ) 2 


2 ^ 1 

'{B/H? S 2 


L(m-i ) 2 (m +|) 2 


_Lu__ _i__ = _ l r~ 


~B/H S nA(BIHy-BIHl%^{BIH) 

For other values of s/B it is necessary to sum the series numerically. The series may be rewritten 

m , „ TO 1 2 


jiff] 


(°’ B’h)- 8 S S E mn F mn 


A r 

f2 S 


[(m +s/B)F m n + (m—$/B)E mn ] n=i n 2 -y/n 2 -\-(s/H) 2 (B/I-jy m 


S [m*-(*/W + ^ a ' (°' 7i’ ll) 

(A3) 


where A r is an arbitrary positive integer and Tin (0 , s/B, B/H) 
denotes a remainder term. In evaluating <r(0, s/B, B/H) in 
the present report, N was usually taken equal to 7 and the 
summations indicated in equation (A3) were carried out ex- 


actly; whereas an approximate value was used for the re- 
mainder term. An approximate formula for 7? v (0, s/B, B/H ) 
is 


n, 




umr(N+l) (A'+i)' [ 


2 


2 


2 -y/ (B/H) 2 (n -b^) 2 -(- 1/4 


B/H ' (71/77) 2 


(b/hy(n+ y ) 


2 ^/(/f///) 2 (l /4)+(-'V + 2 ) J(B/H) 2 +i 


(B/H) 2 (n+£) ’ 0M 2 (A r + 2“) J 

It will now be shown how formula (A4) is obtained. It is easily verified that 

r to''' 

and that 


^(B/H) 2 a 2 +b 2 


(y 2 -f (B/IT) 2 x 2 \ m (BfH)b^ (B/ H) 2 a (B/H) 2 ab 


f" ^£__L 
J 0 7?~2a? 


(A 4) 


(A 5) 
(A6) 


7S1 
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provided a, £>> 0. Now Z? jV ( 0, sjB. U/H) is approximately 
equal to /i iV (0, 0, B/IT) and it is easily found that 


i{ K {o,o,jJiu)=i(±, £ + £ £) 

\» = I in =A r +l- n = A + l m=l/ 


[ n'+mVi/H)*) 3 ' 2 ' 


2 s 




«, (o, |)-4 ( j’-tyj ‘ **+ JJ , 

2 i + 2 + 2 2 

J v+ . dy j x ,i dx ) WhwWF 2+ 

A+ 2 A+ 2 

2 [ 1+ (5/^]J ,~¥ 


Those summations may be approximated by suitable into- If the integrals of equation (A7) are evaluated by means of 


grals so that approximately 


equations (A5) and (Ah), then equation (A4) is obtained. 





APPENDIX B 


LIST OF IMPORTANT SYMBOLS 


B tunnel breadth or diameter 

II tunnel height 

(J tunnel cross-sectional area 

c chord of airfoil or body of revolution 

t maximum thickness of airfoil or body of revolution 

V volume of wing or body of revolution 

S„ maximum' cross-sectional area of body of revolu- 

tion 

x half span of wing 

C n drag coefficient 

S model area on which drag coefficient is based 

U stream velocity 

M Mach number 

B Reynolds immber 

y ratio of specifier heat 'of gas at constan t pressure to 

specific heat at constant volume (c.„/e. r ) 

P mass density 

q dynamic pressure 

A ,K h Ki factors depending on shape of airfoil base profile 
(See equations (11), (12), and (13) and tables 
ff and III.) 

\,Ka,K t factors depending on shape of body of revolution 
(See equations (22), (23), and (24) and tables 
IV and V.) 

h(M') linear compressibility correction factor for virtual 
volume (reference 11) 

t factor depending on tunnel shape and wing-span- 

to-tunnel-breadth ratio (See equations (10) and 
(21) and table I.) 

K total blockage correction (See equation (41).) 

K, c wing-blockage correction (See equations (42) and 

(43).) 

Kb body-blockage correction (See equations (44) and 

(45).) 

K,rk wake-blockage correction (Sec equation (40).) 
Superscript: 

O when pertaining to fluid properties, denotes values 

existing in tunnel far upstream from model; 
when pertaining to airfoil characteristics, denotes 
values in tunnel, coefficients being referred to 
apparent dynamic pressure <]' 


Subscripts: 

(Used only when necessary to avoid ambiguity) 
c denotes values in compressible fluid 

i denotes values in incompressible fluid 

w denotes values for wing 

b denotes values for body of revolution 

wk denotes values for wake. 
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